Emergence of zero-lag synchronization in generic mutually coupled chaotic systems 
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Zero-lag synchronization (ZLS) is acliieved in a very restricted mutually coupled chaotic systems, 
where the delays of the self-coupling and the mutual coupling are identical or fulfil some restricted 
ratios. Using a set of multiple self-feedbacks we demonstrate both analytically and numerically that 
ZLS is achieved for a wide range of mutual delays. It indicates that ZLS can be achieved without 
the knowledge of the mutual distance between the communicating partners and has an important 
implication in the possible use of ZLS in communications networks as well as in the understanding 
of the emergence of such synchronization in the neuronal activities. 
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Two identical chaotic systems starting from almost 
identical initial states, end in completely uncorrelated 
trajectories^, On the other hand, chaotic systems 
which are mutually coupled by some of their internal 
variables often synchronize to a collective dynamical 
behavior 0, The emergence of synchronization plays 
important functioning roles in natural and artificial cou- 
pled systems. One of the most fascinating collective dy- 
namical behavior is the zero-lag synchronization (ZLS), 
known also as an isochronal synchronization. ZLS or 
nearly ZLS was measured in the activity of the brain 
between widely separated cortical regionsjl, 0, 0] , where 
synchronization of neural activity has been shown to un- 
derlie cognitive actsQ. The mechanism of the ZLS phe- 
nomenon has been subject of controversial debate, where 
the main puzzle is how two or more distant dynamical el- 
ements can synchronize at zero-lag even in the presence 
of non-negligible delays in the transfer of information be- 
tween them. 

The phenomenon of ZLS was also experimentally ob- 
served in the synchronization of two mutually chaotic 
semiconductor lasers, where the optical path between the 
lasers is a few orders of magnitude greater than the co- 



herence length of the lasers [9|, [10|, [11|, and the anal- 
ogy between the spiking optical pattern and the neu- 
ronal spiking was also recently established[l3|- This phe- 
nomenon has attracted a lot of attention, mainly because 
of its potential for secure communication over a public 
channel^. In [3] it was recently shown that it is pos- 
sible to use the ZLS phenomenon of two mutually cou- 
pled symmetric chaotic systems for a novel key-exchange 
protocol generated over a public-channel. Note that in 
contrary to a public scheme which is based on mutual 
coupling, private-key secure communication is based on 
a unidirectional coupling [15|, in and it is susceptible to 
an attacker which has identical parameters and is coupled 
to the transmitted signal. The generation of secure com- 
munication over a public channel requires mutual cou- 
pling and was only proven to be secure based on the ZLS 
phenomenon 14 1 . 



alytically that various architectures of coupled chaotic 
maps can exhibit ZLS 17 , • The main disadvantage 
of this phenomenon is that ZLS even between two mu- 
tually coupled chaotic systems can be achieved only for 
very restricted architectures and it is highly sensitive for 
mismatch between the delays of the mutual coupling and 
the self-feedback. These delays have to be identical or 
have to fulfil special ratios. Such a realization might ex- 
ist in a time-independent point-to-point communication, 
but it is far from the realm of communications networks. 

In this letter we first demonstrate the constraint that 
ZLS is achieved only for very restricted ratios between the 
self- feedback and the mutual delays, riTd = mTc, where n 
and m are (small) integers. We next show that one can 
overcome this constraint when multiple self-feedbacks are 
used. For the simplicity of the presentation we mainly 
concentrate on the Bernoulli map, where results of simu- 
lations can be compared to an analytical solution [isl. flit. 
However we observed the reported phenomenon for other 
chaotic maps and systems as well, and it is exempli- 
fied by the ZLS of mutually couple chaotic semiconduc- 
tor lasers, depicted by the Lang-Kobayashi differential 
equations |9ll2d|. 

The cornerstone of our system is the simplest chaotic 
map, the Bernoulh map, f{x) = {ax)modl, which is 
chaotic for a > 1. The dynamical equations of the two 
mutually coupled chaotic units, X and Y, with one self- 
feedback (see solid lines in figure [T]) are given by 



xt = {l- e)f{xt-i) + e[Kf{xt-ra) + (1 - ^^)f{yt-T^] 
yt = {l- e)f{yt-i) + e[«/(yt-rj + (1 - ^),/(yt-rc)] 



(1) 



Recently, it has been shown both numerically and an- 



where and Tc are the delays of the self feedback and the 
mutual coupling, respectively (T^ . The quantities 1 — e, 
SK, and (1 — K)e stand for the strength of the internal 
dynamics, self-feedback and the mutual coupling, respec- 
tively. 

The stationary solution of the relative distance be- 
tween the trajectories of the two mutually coupled 
chaotic Bernoulli maps can be analytically examined 18|, 
[l9l |. Let us denote by 5xt and 5yt small perturbations 
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FIG. 1: A schematic diagram of two mutually coupled units 
at a distance Tc with one self-feedback with a delay equals to 
Td (solid lines). Additional self- feedbacks are denoted by the 
dashed lines. 

from the trajectories xt and yt, respectively. Using the 
ansatz Sxt = c^Sxq and Syt = c^Syo and linearizing equa- 
tions 111, one can find the characteristic polynomial 

c-a{l-e)- aeKc^"^" + ae(l - k^-''-' = (2) 

where A — ln\c\ is the Lyapunov exponent. Simulations 
of the dynamical equations ([1]) and the semi-analytical 
calculation of the maximal Lyapunov exponent of the 
characteristic polynomial (O, indicate that ZLS is the 
stationary solution of the dynamics only when the delays 
of the self-feedback and the mutual coupling fulfil the 
constraint 

nr^j -I- rriTc = (3) 

where the available integers for n,m € Z are functions of 
e and k. Results are exemplified in figure 2 for e = 0.9 
and K — 0.8 (left panel) and for e — 0.9 and k — 0.4 
(right panel). For the left panel, ZLS is achieved for 
the pairs (m, n) = (— l,n) where n = 1,2,. ..,10 and 
(3,— l)[2lj. For the right panel ZLS is achieved for the 
pairs (-1, nU = 1, 4, (3, -1), (5, -1), (7, -1), (3, -2) 
and (5, — 2)[23|. Those lines may have width, so the more 
accurate equation is \nTd — mTc\ < 5, where S » 2. The 




are given by 

N 

xt^{l - e)f{xt-i)+e[K^aif{xt-r^^ )+(l - K)/(j/t_^J] 
1=1 

N 

yHl - £)/(yt-i)+e[K^ai/(2/t-r.,)+(l - «)/(a;t-.J](4) 
1=1 

where N stands for the number of self-feedbacks and the 
parameter ai indicates the weight of the self-feedback 
fulfilling the constraint X^z^i ~ 1- order to reveal 
the interplay between possible which lead to ZLS and 
a given set of {r^j } we first examine in detail the case of 
iV = 2. 

Results of simulations with N = 2 which were con- 
firmed by the calculation of the largest Lyapunov ex- 
ponent obtained from the solution of the characteristic 
polynomial, similar to equation (2), are depicted in fig- 
ure [31 The synchronization points (r^j^ , Td^ ) where ZLS 
is achieved form straight lines. A careful analysis of the 
equations of these lines indicates that their equations are 

niTdi + n2Td^ + TTlTc = (5) 

where ni,n2 and m are integers. The lines may have a 
small width, hence a more accurate equation for the ZLS 
points is \niTd-i + n2Td^ + mrd < 5, where (5 w 2. The 
same equations for the ZLS lines and with similar possi- 
ble width, (5, were obtained in simulations with different 
e, K, ai and Tc, prime and non-prime numbers in the 
range [31,720]. 
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FIG. 2: Simulations and semi-analytic results for the ZLS 
points in the phase space (rd,rc) with a = 1.1, e = 0.9 and 
K = 0.8 left panel and e = 0.9 and k = 0.4 right panel. 

constraint ([3]) indicates that ZLS can be achieved only 
when Tc is accurately known, which is far from the realm 
of communications networks. In order to increase the 
possible ZLS range of Tc for a fixed Td, we added more self- 
feedbacks, as depicted in figure [TJ The generalized dy- 
namical equations for the case of multiple self-feedbacks 



FIG. 3: Simulations and semi-analytic results for the ZLS 
points in the phase space (jd^^Td^ for Tc = 101, a = 1.1,£ = 
0.9, K = 0.8 and cti = 1/2. 

Figure 3 indicates that for e ~ 0.9 and k ~ O.sfl^ , 
for instance, m can take the integers ±1 and ±3 only. In 
order to examine the possible range of the integers \ni\ 
we ran an exhaustive search simulation, —6 < jt.; < 6 and 
TO = ±1, ±3, and obtained integer from equation 
Figure [3] depicts results of such an exhaustive search and 
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the analytical solution of appropriate characteristic poly- 
nomials. The comparison between the results indicates 
the following two main conclusions: (a) n,j ss 6 gives a 
similar synchronization range, (b) the lines have an ex- 
tension of up to 2, hence the actual ZLS points fulfil the 
equation \niTd-^ +n2Td2 +mTc\ < 2 (see the inset of figure 
|4]). Note that a few blue points are missing in the ZLS 
obtained in the semi-analytical solution (red points) in- 
dicating that a few combinations (ni, rt2, m) are missing. 
We also analyze in detail the case of triple self- feedbacks, 
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FIG. 4: ZLS for two mutually coupled Bernoulli maps with 
a — LI, £ = 0.9 and k = 0.8 and with two self-feedbacks, r^i — 
25, Td2 ~ 87. ZLS points obtained from the exhaustive search 
of eq. (5) with m = ±1, ±3 and rii in the range [—6, 6] (blue 
points). The ZLS points obtained from simulation and semi- 
analytical results (red points). The inset is a blow up of a 
section of possible Tc with ZLS. 

equation (4) with iV = 3, and find that ZLS points follow 
the equation \niTd-^ + 712x^2 -|- n^Td^ + rnTc\ < (5 ^ 2, and 
in this case the ZLS points form planes. 

The generalization of the ZLS points for iV = 1, 2 and 
3 to the case of multiple self-feedbacks is 



N 

n-iTdi 

i=l 



(6) 



where and m take bounded integer values. This gener- 
alization was indeed confirmed in simulations and solving 
the characteristic polynomials with up to A'^ = 7. 

In order to obtain a continuous range of Tc for which 
ZLS is achieved, we examined the scenario of 4 different 
Tdi = 11, 15, 18, 150. We select one remarkably large Td 
such that we can see its effect on the range of Tc where 
ZLS is achieved. To measure the quality of the ZLS we 
used the correlation function, which is defined by 



C 



{xtVt) - {xt){vt) 



(7) 



where C = 1 indicates complete ZLS and (...) stands 
for an average over the last 1000 time steps. The cor- 
relation function, C, obtained in simulations is depicted 
in figure 5 and indicates the following results. Multiple 
self- feedbacks result in a continuous range of ZLS for Tc, 



hence it is not required to know exactly the mutual dis- 
tance (value), Tc- Panel (a) of figure 5 indicates that 
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FIG. 5: Simulations results of the correlation, C, as a func- 
tion of Tc for a = 1.1, e = 0.9 and k = 0.8 and four 
Td^ = 11, 15, 18, 150. The weight of the self- feedbacks, ai 
in eq. (4), are in (a) ai = 02 = 03 = 0.25/3 and 04 = 0.75 
and in (b) ai = Q2 = as = 0.65/3 and 04 = 0.35 

there are at least 7 continuous ZLS regimes, each one 
of them is centered at 150n4, where 714 — 0,1, 6 
and the plateaus are extended ~ ±60 around the cen- 
ters (slightly decreases with increasing ^4). This width, 
±60, is much smaller than the ZLS range of the only 
three short self-feedbacks 11,15,18 which was found to 
be ~ [1, 150], indicating that the effective ni,rt2 and 
in eq. (6) are less than 6. This discrepancy is a result 
of the dominated weight of r4 = 150, 04 = 0.75, in fig- 
ure 5(a). For a smaller weight for the largest delay 150, 
Q!4 = 0.35, panel (b) of figure 5, a ZLS is continuously 
achieved up to Tc ~ 700. In this case a weak weight for 
the largest delay results in limited 77.4 which takes the 
values of 0, 1, 2, 3, 4 only, and we expect ZLS in four con- 
tinuous regimes centered around = 0, 150, 300, 450 
and 600 [24] . However these four regimes arc now merged 
by the ±150 width inspired by the strengthened weight 
for the short self- feedbacks, ai = a2 — a-i = 0.65/3[25j. 

In the general case there is an interplay between the 
following three parameters characterizing the set of the 
delay times: Td^^^ which is comparable to Tc, {rd-} <C 
Td^^^ and Aj = Td^^^ - Td^ i = 1,... N - 2, where 
{tcj. } are arranged in an increasing rank order. For 
instance, the following three sets of four self-feedbacks 
(2,6,9,150), (11,15,18,150), (80,84,87,150) are char- 
acterized by the same r^^^^, Ai and A2. What is the 
main difference between the ZLS profile of these sets and 
which set maximizes the continuous range of ZLS? The 
first set opens only a small continuous ZLS regime (~ 20 
for parameters of panel (a)) around 150?t.4, since the time 
delays are very short. The third set almost does not open 
a continuous regime of ZLS, since Td^, Td2, Td^ ^ Ai, A2. 
The maximal continuous ZLS range is achieved when 
short delays r^^ , , r^g are comparable with ~ 6 Ai , 6 A2 
(see eq. (6)) which is a case of the second set. 

Most of the reported simulations were carried out from 
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close initial conditions, however, one can find (e, k) such 
that ZLS is achieved from random initial conditions at a 
comparable time to ZLS with only one time delay, Tc — t^i. 

Similar results were obtained also for mutually coupled 
chaotic logistic maps where the Lyapunov exponent is 
fluctuating in time and is positive only on the average. 

Finally we report that a similar phenomenon of ZLS oc- 
curs in simulations of two mutually coupled semiconduc- 
tor lasers depicted by the Lang-Kobayashi equations [io!]. 
Our simulations are based on the version and the param- 
eters of these equations as in Q, with additional time 
delays. Figure 6(a) depicts the ZLS as a function of 
Tc for the case of 4 time delays Tdi = 3,4,5,20ns and 
Ki — a = ?)Ons~^ and in figure 6(b) for 6 time de- 
lays Trf^ = 11, 12, 13, 14, 15, IQns with ki = a = 257is^^, 
where for both cases the threshold current was p — 1.02. 
For each Tc the duration of the simulation was 7000ns and 
the emergence of ZLS was estimated from the measured 
cross correlation of the last 20 windows of 100nsf26|. Re- 
sults indicate that for the case of 6 delays ZLS is achieved 
in the range ~ [L 80] ns where for the case of 4 time delays 
for ~ [1, 45]ns[27]. These synchronization regimes can be 
explained by equation (6) with n; = 0, ±1,±2 only. It 
is consistent with our simulations of only one time delay 
where ZLS is achieved for = uTc with n — 1,2,3 only 
(instead of 1, 6 for the examined maps) . Note that 
no extension on a time scale of ns is expected, (5 = 0, 
however, preliminary results indicate that a similar phe- 
nomenon occur where — 0.01ns which is comparable 
with the coherence length of the laser. 

The research of LK. is partially supported by the Israel 
Science Foundation. 
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FIG. 6: Simulation results of the correlation, eq. (7), for two 
mutually coupled semiconductor lasers (details in the text). 
Panel (a) for 4 delays = 3, 4, 5, 20ns and panel (b) for 6 
time delays Td, = 11, 12, 13, 14, 15, 16ns. 
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